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A microscopic, stochastic, minimal model for collective and cohesive motion of identical self-propelled particles 
is introduced. Even though the particles interact strictly locally in a very noisy manner, we show that cohesion can 
be maintained, even in the zero-density limit of an arbitrarily large flock in an infinite space. The phase diagram 
spanned by the two main parameters of our model, which encode the tendencies for particles to align and to stay 
together, contains non-moving "gas", "liquid" and "solid" phases separated from their moving counterparts by 
the onset of collective motion. The "gas/liquid" and "liquid/solid" are shown to be first-order phase transitions 
in all cases. In the cohesive phases, we study also the diffusive properties of individuals and their relation to the 
macroscopic motion and to the shape of the flock. 



1. Introduction 



2. Vicsek's model 



The emergence of collective motion of self- 
propelled organisms (bird flocks, fish schools, 
herds, slime molds, bacteria colonies, etc.) is a 
fascinating phenomenon which attracted the at- 
tention of (theoretical) physicists only recently 
f |ll2l4l5l^ l. Particularly intriguing are the sit- 
uations where no "leader" with specific proper- 
ties is present in the group, no mediating field 
helps organizing the collective dynamics (e.g. no 
chemotaxis), and interactions are short-range. In 
this case, even the possibility of collective motion 
may seem surprising. 

However "simple" the involved organisms may 
be, they are still tremendously complex for a 
physicist and his inclination will often be to go 
away from the detailed, intricate, as-faithful-as- 
possible modeling approach usually taken by biol- 
ogists, and to adopt "minimal models" hopefully 
catching the crucial, universal properties which 
may underlie seemingly different situations. 

In this setting, the organisms can be reduced to 
points which move at finite velocity and interact 
with neighbors. This is in fact what Vicsek and 
collaborators did when introducing their minimal 
model for collective motion. 



Vicsek's model |2j consists in pointwise parti- 
cles labeled by i which move synchronously at 
discrete timesteps At by a fixed distance vq along 
a direction 9i. This angle is calculated from the 
current velocities of all particles j within an in- 
teraction range rg, refiecting the only "force" at 
play, a tendency to align with neighboring parti- 
cles: 
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where vl is the velocity vector of magnitude vq 
along direction 6i and ^* is a delta-correlated 
white noise (^ G [— tt, tt]). Fixing ro = 1, At = 1, 
and choosing, without loss of generality, a value 
voAt < ro, Vicsek et al studied the behavior of 
this simple model in the two dimensional param- 
eter space formed by the noise strength 77 and 
p, the particle density. They found, at large p 
and/or small 77, the existence of an ordered phase 
characterized by: 

^^(l(«1)J),>0, (2) 

i.e. a domain of parameter space in which the 
particles move collectively. 

The existence of the ordered phase was later 
proved analytically 7 via a continuous model for 
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Table 1 

Fixed- value parameters used in the simulations. 



the coarse-grained particle velocity and density. 
Vicsek et al devoted most of their effort to study- 
ing the transition to the ordered phase 2 . They 
found numerically a continuous transition charac- 
terized by scaling laws and they tried to estimate 
the corresponding set of critical exponents. 

3. Collective and cohesive motion 

Vicsek's model accounts rather well, at least 
at a qualitative level, for situations where the or- 
ganisms interact at short distances but need not 
stay together. This is for instance the case of 
the bacterial bath recently studied by Wu and 
Libchaber [5]. In this experiment, E. Coli bac- 
teria are swimming freely within a fluid film of 
thickness approximately equal to their size. By 
seeding the system with polystyrene beads and 
recording the trajectories of these passive tracers, 
Wu and Libchaber showed that the bacteria per- 
form supcrdiffusive motion crossing over to nor- 
mal diffusion. We later argued that the supcrdif- 
fusive behavior is likely to be due to the onset of 
collective motion as in Vicsek's model |10lll| . 

When the situation to be described involves the 
overall cohesion of the population, Vicsek's model 
needs to be supplemented by a suitable feature. 
Indeed, an initially cohesive flock of particles will 
disperse in an open space. In other words, no col- 
lective motion is possible in the zero density limit 
of this model. In the following, we extend Vic- 
sek's model to account for the possible cohesion 
of the population of particles. 

The above remark is by no means new. Early 
models for the collective motion of "boids" (con- 
traction of "birdoid", a term used by computer 
animation graphics specialists) do include a two- 
body repulsive- attractive interaction More 
recent works by physicists also included this in- 
gredient, but they either comprised an extra 
global interaction or the actual interaction 
range used extended over the whole flock for the 



sizes considered, making it effectively global 0. 
Another encountered pitfall, from our point of 
view at least, is to enforce the cohesion by the 
confinement to a rather small, close, space [S]. 

Here we want to be, in a sense, in the least- 
favorable circumstances for observing collective 
motion: no leader in the group, strongly noisy 
environment and/or communications, strictly lo- 
cal interactions, and no conflnement at all. 
The "minimal" model presented below is one of 
the simplest possible ones satisfying these con- 
straints. 

4. A minimal model 

In addition to the possibility of achieving cohe- 
sion, we also want to confer a "physical" extent 
to the particles, a feature absent from Vicsek's 
point-particles approach. Adding a Lennard- 
Jones-type body force / acting between each pair 
of particles within distance tq from each other 
offers such a possibility. 

Equation (Q) is then replaced by 



gt+l 



= arg 



(3) 



where a and (3 control the relative importance of 
the two "forces" . The precise form of the depen- 
dence of the body force on the distance between 
the two particles involved is not important. It 
is enough to ensure a hard-core repulsion at dis- 
tance Tc and an "equilibrium" preferred distance 

In the following, we use 



fij — e 




if Tij < rc , 

if Tc < nj < Ta 

ff 7-a < nj < To 



(4) 



with rij the distance between boids i and j, Cij 
the unit vector along the segment going from i to 
j, and the numerical values Tc = 0.2, = 0.5 and 
ra = 0.8. 

We have also tested other types of noise term in 
the model. In particular, considering the noise as 
the uncertainty with which each boid "evaluates" 
the force exerted on itself by the Ni neighboring 
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Figure 1. Cohesive flocks of 128 particles in a square box of linear size 32 with periodic boundary 
conditions (for parameters see TableCJ. (a): immobile "solid" at a = 1.0 and /3 = 100.0 (20 timesteps 
superimposed), (b): 3 snapshots, separated by 120 timesteps, of a "flying crystal" at a = 3.0 and 
f3 = 100.0. (c): fluid droplet (a = 1.0, (3 — 2.0, 20 consecutive timesteps). (d): moving droplet (a = 3.0, 
f3 = 3.0, 20 consecutive timesteps). In (b) and (d), the arrow indicates the (instantaneous) direction of 
motion. 
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Figure 2. Order parameters at p = 1/16, L — 128, a = 1.0. (a) : "gas/liquid" transition, inset: cluster 
mass distribution at the coexistence point [3 — 1.0. (b) : "liquid/solid" transition, inset: pdf of order 
parameter A for the "hquid/sohd" transition, /3 = 36 plain line ("liquid" phase), /3 = 48 dotted line and 
/3 = 60 dashed line ("solid" phase), (other parameters as in Table QJ. 



boids leads to change Eq. lO to: 



qt+l 



arg 



(5) 



where u* is a unit vector of random orientation. 
There are delicate issues related to the choice of 
the noise term, in particular with respect to the 
critical properties of the transition to collective 
motion 12 . We mostly considered, in the fol- 
lowing, the noise term as prescripted above in 
Eq. 0. A detailed analysis of the influence of 
the nature of the noise is left for future studies, 
but we are confident that the results presented 
here hold generally, at least at a qualitative level. 

Finally, in order to ensure that each parti- 
cle only interacts with its "first layer" of neigh- 
bors (within distance j'o), we calculate, at each 
timestep, the Voronoi tessellation of the popu- 
lation — this has the additional advantage of 
providing a natural definition of the "cells" asso- 
ciated with each particle. The interacting neigh- 
bors are then retricted to be those of the particles 



within distance ro which are also neighbors in the 
Voronoi sense. 

5. Typical phases 

One can easily guess the "phases" that the 
above model can exhibit for a fixed noise strength 
77 (77 = 1.0 in the following, for a summary of pa- 
rameters see Table ^ . We now present them in 
a qualitative manner. The results presented be- 
low were all obtained in the two-dimensional case, 
but most of them hold in three dimensions. When 
the body force is weak (small (3 values), the cohe- 
sion of a fiock cannot be maintained. In a finite 
box (finite particle density p), one is left with a 
gas- like phase (not shown). An arbitrarily large 
flock in an infinite space disintegrates, eventually 
leaving isolated random- walking particles. 

For large enough /3, we can expect the co- 
hesion to be maintained. Figure ^ shows such 
cohesive flocks. The internal structure of the 
flocks depends also on (3: for large body force, 
positional quasi-order is present, and the parti- 
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Figure 3. Mean square distance between initially- 
neighboring particles vs time for a flock of 10000 
boids {L = 400, p = 1/16) in logarithmic scales, 
(a) non-moving cohesive droplet {a — 1.0, (3 — 
25.0, 35.0, 45.0 and 150.0 from top to bottom, the 
dashed line has slope 1); (b) moving cohesive flock 
(a = 3.0, P = 40.0, 55.0, 75.0 and 150.0 from top 
to bottom). Asymptotic transition points were 
measured at a = 1.0, /3ls = 45.3 and at a = 3.0, 
/3ls — 76.1 (see below). Other parameters as in 
Table ^ Insets: data collapse from which the 
waiting time r can be estimated. 



cles locally form an hexagonal crystal (Fig. QJ^b) . 
For intermediate values of f3, no positional order 
arises, and the flock behaves like a liquid droplet 
(Fig.DJd). 

The influence of the alignment "force" is mani- 
fested by the global motion of the flock: for large 
enough a values, the flock moves {V > 0). De- 
pending on /3, one has then either a "moving 
droplet" (Fig.[TJi) or a "flying crystal" (Fig.IB)). 

6. Order parameters 

Order parameters have to be defined to allow 
for a quantitative distinction between the phases 
described above. 

The limit of cohesion separating the "liquid" 
phases from the "gas" can be determined by 
measuring the distribution of the sizes of par- 
ticles clusters, thanks to an implementation of 
the Hoshen-Kopelman 15 algorithm. A cohe- 
sive flock is then one for which n, the size of the 
largest cluster, is of order N, the total number of 
particles. Below, we use the criterion n/N = ^ to 
define the transition. Increasing (3, n/N sharply 
rises to order-one values (Fig.[21i). 

The "liquid/solid" transition takes place when 
P is large enough so that cohesion of the pop- 
ulation is ensured. To determine this onset of 
positional quasi-order within a finite but arbi- 
trarily large cohesive flock, we first observe that, 
whether in the collective motion region or not, the 
"liquid" and "solid" phases can be distinguished 
by the fact that particles diffuse with respect to 
each other in the "liquid" , whereas neighboring 
particles always remain close to each other in 
the "solid" (Fig. E)). To be more precise, in the 
"liquid" , initially close-by particles remain so for 
some trapping time r (which can be defined or 
used to collapse the curves of Fig. ^ . Approach- 
ing the "solid" phase (by increasing /3), r di- 
verges. In addition, since we are not dealing with 
a translation-invariant system, we should also 
distinguish between the diffusion properties of 
the particles depending on their relative position 
within the flock. We have thus defined different 
"sectors" ("core", "head", "tail", "sides") as ex- 
plained in Fig. ^ If all sectors are roughly equiv- 
alent in non-moving droplets (at least sufficiently 



6 




-20 -10 10 20 



Figure 4. Snapshot of a moving flock of 4096 
boids (coordinates centered on the position of the 
center of mass (CoM)). The arrow indicates the 
instantaneous direction of motion of the flock. 
The solid circle, centered on the CoM, has a ra- 
dius equal to the root-mean-square of all boids' 
distances to the CoM. The circle with a 4 times 
smaller radius deflnes the "core" region (flUed cir- 
cles near the CoM) . The "head" contains all boids 
outside the larger circle which are also within a 
cone of opening angle 7r/4 centered on the direc- 
tion of motion. The "tail" region is deflned in an 
opposite manner. 



far from the "liquid/solid" transition), some dif- 
ferences are observed within moving droplets: the 
outer regions and in particular the head are more 
active, whereas cohesion is stronger in the core 
(Fig. 0). Thus, in principle, different trapping 
times r can be deflned for the different regions 
of a (moving) flock. But the depth of the outer, 
"more liquid" , layer does not depend on the flock 
size (if it is big enough), so that, in large flocks, 
most of the population behaves as the core. Con- 
sequently, the relative diffusion averaged over the 
whole flock suffers from finite-size effects, but 
they disappear in the large-size limit (see below). 
To sum up, the trapping time r (measured on 
all particles of the flock) is a good quantity to 
track the "liquid/solid" transition. However, in- 
stead of directly estimating of r, we measured A, 
the relative diffusion over some large time T of 
initially-neighboring particles: 




where the ni particles j are the neighbors of parti- 
cle i at time t. Time T is taken to be proportional 
to the volume of the system, which ensures that 
A records, in the large-size limit, an asymptotic 
property of the system (because T :s> r) . Clearly, 
A ~ 1 in the "liquid" phase for a large enough 
system, while A in the "solid" phase. In- 
deed, increasing [3, A falls off sharply (Fig. [^b). 
The transition point is chosen to be at A = i. 

Finally, we must be able to distinguish the 
regimes for which collective motion arises. To 
this aim, we use the average velocity V , as de- 
fined in Increasing a, V reaches order-one 
values. The onset of collective motion is chosen 
to be at F = ^0/2. In the gas phase, one ex- 
pects V — Q independently of the strengh of the 
alignment force. Nevertheless, this phase is rather 
sharply divided in two when studying the model 
at finite particle density. The largest cluster size 
may then be small {n <^ N), but it is (almost) 
always larger than one. At any given time, thus, 
n > 1, and the collective motion order parame- 
ter V can be defined restricted to the particles 
belonging to the largest cluster. (Since clusters 
merge and break, the particles involved generally 
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Figure 5. Growth of the mean square distance 
between initially-neighboring particles depending 
on their position within a cohesive "liquid" flock 
of 10000 boids (p = 1/16, /3 = 40). (a) a = 1.0, 
cohesive non-moving droplet; (b) a = 3.0, moving 
droplet. Solid line: all boids in flock. Dotted 
line: core region, dashed line: head region, dash- 
dotted line: tail region (as defined in Fig.^J) Insets 
: diffusion within flock head in linear scales, the 
solid line is just a guide to the eye. 



change along time.) 

7. Phase diagram 

After a brief discussion of the nature of the 
transitions involved, we first present the phase 
diagram of our model for a finite density of par- 
ticles in a large and fixed box size. Then we esti- 
mate finite-size effects on the location of the phase 
boundaries. Finally, we argue that the phase dia- 
gram can also be defined in the zero-density limit 
where an arbitrarily large flock wanders in an in- 
finite space. 

7.1. Nature of the transitions 

As expected in usual phase transitions, we 
found that in our model the "gas/liquid" and "liq- 
uid/solid" transitions are first-order. In insets of 
Fig.Et and b, we show the evolution of the prob- 
ability distribution function of the order param- 
eter as one crosses these transition lines. The 
bimodal character of these pdf at the transition 
is typical of first-order phase transitions, indicat- 
ing the coexistence of two metastable states. At 
the "gas/liquid" transition point, dispersed and 
aggregated boids coexist and there are exchanges 
between the two phases along time. At the "liq- 
uid/sohd" transition point, cohesion is ensured 
and one observes the quasi-frozen regions in an 
otherwise more "liquid" flock. The quasi-solid 
parts are often located in the core of the flock. 
(Fig.©. 

The nature of the transition for the onset of 
collective motion is a delicate issue in our model. 
Whereas its second-order character is rather well- 
established for Vicsek's core model, we recently 
discovered that, in fact, the implementation of 
the noise term may change the nature of the tran- 
sition. The detailed investigation of this, in par- 
ticular in presence of the cohesive force, will be 
presented in a future publication |12| . 

7.2. Fixed population and fixed box size 

A systematic scan of the (a, j3) parameter plane 
was performed for a flock of iV = 2025 particles 
living on a square surface of linear size L = 180 
{p = 1/16) with periodic boundary conditions. 
Using the criteria defined above, we obtained the 
phase diagram presented in Fig. 
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Figure 6. Short-time trajectories of freezing 
droplets of = 512 boids (5000 timesteps are 
shown, the motion of the center of mass and the 
sohd rotation around it have been substracted) . 
(a): in the moving phase a = 3.0, /3 = 70: the 
inner part appears more sohd, while the head is 
clearly more liquid (the arrow indicates the in- 
stantaneous direction of motion) . (b) : in the non- 
moving phase a = 1.0, /3 = 50: one can distin- 
guish an outer liquid layer from the almost solid 
core. 



For each parameter value, we used an initially- 
aggregated flock. We let the system evolve during 
a time t a U^, and then we recorded each or- 
der parameter and its histogram along time. The 
transition points were determined by dichotomy, 
the precision of which is reflected in the error 
bars. 

The basic expected features are found: the 
horizontal "gas/liquid" and "liquid/solid" tran- 
sitions arc crossed by the vertical "moving/non- 
moving" line. Near this line, however, one ob- 
serves a strong deformation of the "gas/liquid" 
and "liquid/solid" boundaries. This cannot be 
understood without a careful study of the collec- 
tive motion transition 12 . 

Note also that the "gas" phase itself is crossed 
by the line marking the onset of collective motion, 
using, as explained above, the average velocity V 
of the n particles of the largest cluster as the order 
parameter. 

7.3. Finite size and saturated vapour ef- 
fects 

We are ultimately interested in the possi- 
bility of collective and cohesive motion for an 
arbitrarily-large flock in an infinite space. The 
phase diagram of Fig. |3 was obtained at a fixed 
system size and constant density. Thus both lim- 
its of infinite-size and zero-density have to be 
taken to reach the asymptotic regime of interest. 
Of course this is mostly relevant to the onset of 
cohesion (the "gas/liquid" transition). Here we 
first study each limit separately, i.e. we inves- 
tigate finite-size effects at fixed particle density 
and expansion at fixed particle number. Then we 
discuss the double- limit regime of interest. 

Performing such a task for the whole parameter 
plane far exceeds our available computer power. 
We restricted ourselves to three typical cases: in 
the non- moving phase {a — 1.0), in the moving 
phase {a = 3.0) and near the transition to collec- 
tive motion {a ~ 1.75 for the "gas/liquid" tran- 
sition and a = 2.1 for the "liquid/solid" transi- 
tion). 

In simulations performed at p — 1/16, vary- 
ing N and L, the transition points Pqi^{N) and 
P^c,{N) converge exponentially as L ~ \/N/p 
increases. In Fig. |S| we show this for both the 
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Figure 7. Phase diagram at p = 1/16, i = 180 
(other parameters as in Table QJ. S : sohd, MS : 
moving solid, L : liquid, ML : moving liquid, G : 
gas, MG : moving gas. Dashed line : transition 
line of collective motion. 



"gas/liquid" (a) and the "liquid/solid" (b) tran- 
sitions for a = 3 (moving phases). This al- 
lows to determine asymptotic transition points. 
Note that exponentially-decreasing finite-size ef- 
fects are typical of first order phase transition at 
equilibrium |16| . In the two other cases (non- 
moving phase, and onset of motion), the results 
are similar, with the asymptotic values being 
quickly reached in the non- moving phase (a ~\) 

m 

In simulations performed at fixed N varying 
L, we study instead the expansion of the system. 
The (finite) spatial extent induces a confinement 
effect which increases the pressure at the coexis- 
tence point. We thus expect a displacement of the 
transition point. We find that /3gl(p) also con- 
verges exponentially, and thus transition points 
/3gl(^) are well-defined. In Fig.lHt, we show this 
only for the "gas/liquid" transition at the a val- 
ues used in Fig.|S|i, since the "liquid/solid" tran- 
sitions points were observed to be independent of 
the box size. 

7.4. Zero-density limit 

The above results provide evidence that our 
model possesses well-defined, asymptotic phase 
diagrams at either fixed particle density or at 
fixed number of particles. The double limit men- 
tioned above can be approached in essentially 
three different ways. 

One can take one limit after the other one, re- 
peating either the calculations of Fig. at lower 
and lower densities, or those of Fig. |SJi at larger 
and larger flock size. However, this straightfor- 
ward program involves very heavy numerical sim- 
ulations. Therefore, we only considered two cases 
(a = 1.75 and a = 3.0). As expected, the tran- 
sition points converge (exponentially) indepen- 
dently of the order with which the two limits are 
taken, yielding estimates of the zero-density limit 
(Fig. inj. These estimates are compatible with 
each other. There are three different sources of 
error: statistical and systematic errors when eval- 
uating the order parameters of each system, and 
then fitting errors when determinating /3gl(-^) 
and /3gl(p)- For a= 1.75, we find a 9% difference 
in the estimated asymptotic values whose origin is 
probably the statistical errors on the larger flocks 




Figure 8. Finite size effects on transitions at a fixed density or at fixed number of particles at a = 3.0 
(moving phases). The dashed lines are the asymptotic values determined by exponential fits (solid 
lines), (a) variation of the "gas/liquid" transition point /3gl at p = 1/16 us linear size, (b) variation 
of the "liquid/solid" transition point /3ls at p = 1/16 vs linear size, (c) /3gl vs inverse density 1/p for 
N = 4096. 



simulations. 

7.5. Evaporation of a flock 

There exists a third manner of approaching the 
zero-density limit of the "gas/liquid" transition. 
It consists in quenching a cohesive flock observed 
at large-enough /3 to a lower /3 value. If the flock 
is quenched below the "gas/liquid" line, it will 
"evaporate" . The largest cluster will progres- 
sively lose particles before finally equilibrating in 
the gas phase. This transient can be expected to 
be governed by an effective surface tension 7 and 
the boundary of the largest cluster should then be 
governed by some AUen-Cahn law |17) : u„ = 7K 
where w„ is the (local) normal velocity and k the 
local curvature. Assuming that the mass and the 
surface of the main cluster remain proportional, 
in the two-dimensional problem: n — piocTri?^, 
we can integrate this equation, and we obtain the 
relation : 

n = N — Xt with A = 7rpioc7 

where pioc is the local density: the size of a circu- 
lar flock should decrease linearly in time, with the 
proportionality constant providing an estimate of 



A. 

This is indeed what can be observed in our 
model fFig.llOh, and b). In these experiments, an 
initially large cohesive flock is prepared at some 
(3 value above the transition. The system size 
is taken so as be close to the zero-density limit 
(p = 1/256 in Fig. We measured Tov, the 

time taken by the largest cluster to reach a given 
normalized mass n/N, as well as its surface at the 
same mass/time. We thus checked the propor- 
tionality between mass and surface and between 
mass and time, after transients and before the 
system approaches the equilibrium (Fig. llOb and 
b). Note that due to the abrupt change of pa- 
rameters the flock actually first expands after the 
quench before setting in the "true" evaporation 
regime (Fis. llOb). This experiment can, at first 
sight, be thought of being free of confinement ef- 
fects, and the AUen-Cahn law is expected to be 
satisfied at all times. 

Figure [THfe shows the dependence of X on (3: X 
quickly decays and reaches very small values for 
/3 ~ 0.8. On a logarithmic scale (Fig. we 
can distinguish two exponential regimes, on each 
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Figure 9. Zero-density limit of the cohesion tran- 
sition for (a) a = 1.75 and (b) a = 3.00 (other 
parameters as in Table Transition points at 
zero-density for different system sizes N (circles), 
and at infinite-size for different (inverse) densities 
(squares). The solid lines are exponential fits. 
At a = 1.75, I3gl{N) 11.7 and Pgl{p) ^ 12.7 
; at a = 3.0, PohiNorp) 5.6. 



side of this value, which corresponds roughly to 
the finite-size threshold Pqi^{L) determined above 
(and is rather far from the asymptotic threshold 
determined above to be around 1.5) A simple ar- 
gument can account for this behavior: Consider 
a flock of N boids and suppose there is no short- 
time expansion (such as seen in Fig. llOb ^. The 
mass of the largest cluster n(t), from then on, de- 
creases linearly with time until it reaches, at time 
Too, the equilibrium value neq. At every time t, 
we have: 

N~n{t)^^{N-n,^) = Xt. (7) 

From the theory of first order phase transitions 
of systems at equilibrium (see for instance |lf)|'l. 
we expect the order parameter Ucq/N to behave 
like 

^^i[l + tanh(XiL'^(/3-/3GL))] , 

in the vicinity of the transition, where Ki is a 
constant. Moreover, we expect that Too depends 
linearly on N (see Fig. llOb . where results for two 
different system sizes have been superimposed). 

From a mean-field point of view, Too/iV can be 
interpreted as the mean time required for a parti- 
cle to escape from the interaction of another bold. 
Given the interaction we use (see Eq. |0J), we can 
assume that the potential is harmonic, so that the 
escape time is proportional to the exponential of 
the potential depth (3. Finally, we get 

Too V N J ' 
A cc exp (-i^2/3) [1 - tanh [K^L^ {fi - /3GL(i))]] 

Approximating tanh by (— 1 + exp) and (1 — exp) 
sufficiently far below and above /3 — PchiL), 
we find the two exponential regimes mentioned 
above. More precisely, for f3 < /?GL(i), the sur- 
face tension should be almost independent on sys- 
tem size, whereas for f3 > PchiL), A is governed 
by the finite size effects and its slope increases 
like L'*. Quantitative agreements with the above 
approximation would require too large numeri- 
cal calculations, but our partial data is consistent 
with the above predictions. 
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Figure 10. Evaporation at p = 1/256, a ~ 1.0 (other parameters as in TableCJ. (a) : time vs normalized 
cluster mass, L — 512, and 1024, /3 — 0.2. (b) : surface vs normalized cluster mass, L — 512, /3 — 0.2. (c) 
surface tension vs /3, in linear scale, (d) same as (c) in log-lin. scales. The dashed lines are only a guide 
to the eye. 
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Figure 11. Mean-square displacement of the cen- 
ter of mass vs time for a cohesive droplet (solid 
lines) and for a crystal (dashed lines), N — 32, 
p = 1/64, logaritimic scales, (a): in the non- 
moving phase, a = 1, /3 = 40 ("liquid") and 
65 ("solid"), (b): in the moving phase, a = 3, 
/3 = 40 ("hquid") and 95 ("solid"). The straight 
lines have slope 1 or 2. 



A concluding remark to this investigation of 
the effective surface tension governing the evapo- 
ration of a flock is that, contrary to naive argu- 
ments, this method does not offer much advan- 
tage over the double-limit procedure presented 
in Section 17.41 Indeed, as shown above, it only 
allows a rather easy determination of Pgl{L), 
while the finite-size effects remain hard to esti- 
mate quantitatively. 

8. Micro vs macro motion 

The existence of cohesive phases being now 
well-established, a natural question is that of the 
properties of the trajectories of cohesive flocks 
(the "macroscopic" motion) and it is interest- 
ing to compare those to the trajectories of the 
individuals composing the flock ("microscopic" 
motion). Postponing again the account of what 
happens in this respect near the onset of col- 
lective motion to a further publication we 
studied, for the four possible cohesive phases, the 
mean square displacement (Ar'^it)) of the center 
of mass, as well as of individual boids. 

Our model being essentially stochastic, it is 
no surprise that, at large times, we observe that 
(Ar^(t)) ^ t, i.e. the flock performs Brownian 
motion, in all cases fFig. 1111 and I12|l . When in 
a moving phase (either "liquid" or "solid"), this 
random walk may consist of ballistic flights sepa- 
rated by less coherent intervals during which the 
flock often changes direction. This is testified by 
the baUistic part ((Ar^(t)) ~ i^) of the plots in 
Fig. II lb and by the trajectories itselvesEl The 
moving and non-moving phases can also be dis- 
tinguished by opposite finite-size effects: in the 
non- moving phases, the diffusion constant of the 
macroscopic random walk decreases with system 
size (like l/N), whereas the ballistic flights' du- 
ration increases with system size in the moving 
phases (Fig. 

Comparing the macroscopic motion of "liquid" 
and "solid" flocks is not well defined since, in the 
phase diagram, the line marking the onset of col- 
lective motion is not straight. Nevertheless, at 
comparable distances from this line, one notices 
that ballistic flights tend to be longer for flying 
crystals than for moving droplets. Similarly, the 




Figure 12. Trajectories of the center of mass in the non-moving phase ((a) : a = 1) and in the moving 
phase ((b) : a = 3). Flock of A'' = 32 boids, p = 1/256, (3 = 40.0. At long times (a,b) 10^ timesteps. At 
short times (insets): 1000 timesteps shown. 



diffusion constant of "solid" , non-moving flocks is 
smaller than that of non-moving droplets. 

Finally, even though we have already con- 
sidered the mutual dispersion of initially- 
neighboring boids (see Fig. O, we also studied 
the diffusive properties of individual boids within 
cohesive flocks, substracting out the translation 
motion of the center of mass. In "solid" flocks, 
one can hardly record any such motion. For 
droplets (Fig. I14|l . on the other hand, this micro- 
scopic motion depends on the macroscopic mo- 
tion. When the droplet is fixed, the diffusion 
is normal (Ar^) ^ t, whereas a bold which be- 
longs to moving flock diffuses as (Ar^) ^ i", with 
a ^ 4/3. We interpret this as being due to some 
mesoscopic "hydrodynamical" structures within 
moving flocks (jets, vortices, etc.). J. Toner and 
Y. Tu [H] have shown via a mesoscopic equation 
that collective motion induces transverse correla- 
tions even in a co-moving frame. Therefore, bold 
diffusion must be faster than Brownian. They 
predicted an exponent equal to 4/3, with which 



our results are in good agreement (Fig. 1141 top 
line). 

9. Summary and Perspectives 

We have introduced a simple model for the col- 
lective and cohesive motion of self-propelled par- 
ticles. We have described its various dynamical 
phases, defined order parameters to distinguish 
them, and presented a typical phase diagram at 
large but finite number of particles N and large 
but finite system size L (Fig. [7)). Even though we 
have provided evidence that this phase diagram 
possesses well-defined N ^ oo and L —t oo limits, 
these limit diagrams require too heavy numerical 
simulations to be determined at this stage. 

We have also argued that the double limit of 
an arbitrarily large flock evolving in an infinite 
space is also well-defined. Here also, the map- 
ping of the phase diagram in this limit is cur- 
rently out of reach. But the existence of cohe- 
sive phases in a model of noisy, short-range inter- 
action, identical particles is ensured, which was 
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Figure 13. Mean square displacement of the cen- 
ter of mass vs time for a moving droplet of size 
N = 1024, 256, 128, and 32 from top to bottom 
(logarithmic scales, p ~ 1/64, a — 1.95, (i — 40). 
The solid lines have slope 2 and 1 . The larger the 
flock the more ballistic its motion at short times. 
The inset represents the cross-over time between 
the ballistic and the brownian motion. 
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Figure 14. Internal mean-square displacement vs 
time of the individual boids of a cohesive droplet 
of 10000 boids (p = 1/16, /3 = 40, the motion of 
the center of mass has been substracted) . Lower 
curve: non-moving case (a = 1), normal diffu- 
sion (the associated dashed line has slope 1). Top 
curve: moving case (a — 3), superdiffusion (the 
associated dashed line has slope 4/3, as predicted 
by Toner and Tu 0). 



one of our primary goals. Futhermore, we can 
sketch the zero-density asymptotic diagram from 
our partial knowledge (Fig. A few remarks 
are in order to explain the expected shape of the 
asymptotic "gas/liquid" boundary: in the non- 
moving phase, this transition is almost indepen- 
dent on a and size-effects are negligible. We thus 
expect this line to be horizontal, in agreement 
with mean-field arguments (see Appendix). Sim- 
ilarly, our preliminary study of the onset of collec- 
tive motion shows that one can go directly from 
a cohesive non-moving droplet to the incohesive 
phase and that in this region the onset of collec- 
tive motion is roughly independent on /3, yielding 
a vertical boundary. 

More work should be devoted to the determina- 
tion of the asymptotic phase diagrams mentioned 
above, as well as to a quantitative study of the 
onset of collective motion. Even in the simplest 
case of the non-cohesive Vicsek-type models, it 
can be second or first order, depending on the 
nature of the microscopic noise in the model [H] . 



In both cases, we have started to uncover a rich 
interplay between collective motion, critical fluc- 
tuations, rotation modes, and shape dynamics in 
the transition region. 

Not even mentioning the study of our model 
in three dimensions and its possible applications 
to particular real-world situations (e.g. biology, 
zoology, and robotics) the work presented here is 
probably only the beginning of the exploration of 
this new type of non-equilibrium systems. 

Appendix: Mean-field approach 

To understand the interplay between alignment 
and body force, we simplify our system via a 
Hamiltonian model. The total energy is the 
sum of kinetic energy, two-body interaction en- 
ergy, and a term related to the velocity alignment 
"force" . We thus write 
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where we have introduced an external field hg. 

In the spirit of the mean-field approach, we as- 
sume that the fluctuations are negligible so that 
each component can be integrated independently 
and approximated by its averaged value. Thus 
the contibution from the two-body interaction be- 
comes: 

1 



U (ri, ...,rAr) 



2 ^ 



u{\n -rj\) 



Un = 



I3pu{r)(fr 2apl3, 



where a is a constant which depends on the po- 
tential. The alignment energy is computed with 
the coarse-grained velocity = (u) and the aver- 
ts : 



age number of neighbors pirr^ 

{ho- 



N 



N 



apsif 



N 



In the following, we use /iq + oips(p = h. 

Because of the separation of phase space vari- 
ables, each term of the Hamiltonian contributes 
a factor to the partition function: 
• the two-body coupling : 



{S-hNf e^v[aNp-^ 



(8) 



where b = ^i?^ is the hard-core surface and S is 

the whole surface, 

• the alignment coupling : 

2nmv,lJ-^\Y , (9) 



where Iq is a Bessel function. The free energy of 
the system is thus written: 



F 



NkT ^ ^\ h ) 2kT 



1 + ln {p-^ -b) 



ap 



kT 



Inio 



kT 



psif 



2kT 



Imposing to be at a minimum of the free energy, 
we find that the average velocity is the solution 





Figure 15. Sketch of the asymptotic phase dia- 
gram in the zero density limit. Filled circles in- 
dicate points determined numerically as in Sec- 
tion!^] See text for more details. 
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of a self-consistent equation: 



f_ 



Ii / voh 
lo \kT 



and that the critical point is defined by 
svQpa — 2kT . 



(10) 



(11) 



Collective motion thus emerges via a continous 
phase transition whose critical exponents are 1/2 
for the order parameter and 1 for its susceptibil- 
ity, independently of the cohesive force. 

Furthermore, the necessary concavity of the 
free energy function implies a first order cohe- 
sion transition, since the second derivative of the 
free energy has two zeros. Their positions define 
the phase coexistence region. We determined the 
stability limit of the gas phase. As positions and 
velocities are decorrelated, we studied two cases : 
first without any motion, and then with collective 
motion. In the first case: 



kT 



2ap{l - bp) 



(12) 



which means that the transition line does not 
depend on the alignment parameter (Fig. I16() . 
Note also that there is no transition point in the 
zero-density limit of this model. We numerically 
solved the case with collective motion and found a 
stabilization of the "liquid" phase (Fig. IK) |) . This 
is an effect of the assumption that velocity fluc- 
tuations are negligible. At the onset of motion, 
we expect that fluctuations of the macroscopic 
velocity diverge, leading to strong density fluctu- 
ations. This could explain the de-stabilization of 
the "liquid" phase (Fig. I15() and its stabilization 
in the mean field model 1161 where such fluctua- 
tions are by dcflnition absent. 
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